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ABSTRACT. The paper deals with the semi-classical behaviour of quantum dy- 
namics for a semi-classical completely integrable system with two degrees of free- 
dom near Liouville regular torus. The phenomomenon of wave packet revivals is 
demonstrated in this article. The framework of this paper is semi-classical analy- 
sis (limit : h — > 0). For the proofs we use standard tools of real analysis, Fourier 
analysis and basic analytic number theory. 



1. Introduction 

1.1. Motivation. In quantum physics, on a Riemannian manifold (M,g) the evo- 
lution of an initial state ipo £ L 2 (M) is given by the famous Schrodinger equation 



= PhW); 0(0) = 

Here/z > is the semi-classical parameter and the operator Pj, : D (P/,) C L 2 (M) — > 
L 2 (M) is /z-pseudo-differential operator (for example P/, = — + V). In the 
case of dimension 1 or for completely integrable systems, we can describe the 
semi-classical eigenvalues of the Hamiltonian P/, and by linearity we can write 
the solutions of the Schrodinger equation. Nevertheless, the behaviour of the so- 
lutions when the times t evolves in larges times scales remains quite mysterious. 

In dimension 1, the dynamics in the regular case and for elliptic non-degenerate 
singularity have been the subject of many research in physics [Av-Pe], [LAS], 
[Robil], [Robi2], [BKP], [Bl-Ko] and, more recently in mathematics [Co-Ro], [Rob], 
[Paul], [Pau2], [Lab2]. The strategy to understand the long times behaviour of dy- 
namics is to use the spectrum of the operator P/, . In the regular case, the spectrum 
of Ph is given by the famous Bohr-Sommerfeld rules (see for example [He-Ro], 
[Ch-VuN], [Col]) : in first approximation, the spectrum of P/, in a compact set is a 
sequence of real numbers with a gap of size h. The classical trajectories are periodic 
and supported on elliptic curves. Always in dimension 1, in the case of hyperbolic 
singularity we have a non-periodic trajectory. The spectrum near this singularity 
is more complicated than in the regular case. In [Lab3] we have an explicit de- 
scription of the spectrum for an one-dimesional pesudo-differential operator near 
a hyperbolic non-degenerate singularity. The article [Lab4] deals with the quan- 
tum dynamics for the hyperbolic case. So, in dimension 1, we get the full and 
fractionnals revivals phenomenon (see [Av-Pe], [LAS], [Robil], [Robi2], [BKP], 
[Bl-Ko], [Co-Ro], [Rob], [Paul] for the elliptic case and see [Lab4] or [Pau2] for 
the the hyperbolic case). For an initial wave packets localized in energy, the dy- 
namics follows the classical motion during short time, and, for large time, a new 
period T rev for the quantum dynamics appears : the initial wave packets form 
again at t = T rev . 
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Physicists R. Blhum, A. Kostelecky and B. Tudose are interested in the case of 
the dimension 2 (see [BKT]). Our paper presents some accurate results on the time 
evolution for a generical semi-classical completely integrable system of dimension 
2 with mathematical proofs. 



1.2. Results and paper organization. Here the quantum Hamiltonian is of the 
type P/, = F (Pi, P2) where F is a real polynomial of two variables and Pi, Pi are 
semi-classical one dimensional harmonic oscillators (see section 2 for details). By 
a diffeomorphism this Hamiltonian is less particular than it seems to be, since 
it gives the spectrum of any completely integrable system with two degrees of 
freedom near regular torus or around elliptic singularity [VuN]. Therefore, the 
Hamiltonian study leads to a study more or less general but which is not obvious 
in dimension 2. In this paper, we consider an initial state ipo localized near some 
regular Liouville torus of energies (E\, E2) and we study the associated quantum 
dynamics. To understand the behaviour of dynamics, we interested in the evolu- 
tion of the autocorrelation : 

a (0 = |fa(t)/Vto)tfOR2) 

Due to the simple nature of the Hamiltonian operator the autocorrelation function 
can be write as a serie : 



(0 



-\-OQ -\-OQ 

E E K m |V^(w,«) 



n=0 m=0 



where t„ = <x>\h \ n + |j , ]i m = <^2^ (jn + 5") are eigenvalues of the one-dimensionnal 
harmonic oscillators Pi, P?. The sequence (a n ,m) n m is just the decomposition of the 
initial vector ipo on the Hermitte's eigenbasis of L 2 (]R 2 ). 

Most of the paper (section 3 and 4) consists in estimating and analyzing the 
function a(f) for large times scales (f < l/h s with various s > 0). We use Taylor's 
formula to expand the phase term tF (r n , ]i m ) /h in the variables (n,m); first in 
linear order (section 3), then to quadratic order (section 4). 

In the section 3, we study the linear approximation &\(t) (see definition 3.6) of 
the autocorrelation function, valid up on a time scale [0,1 /h a ] where < a < 1. 
The dynamics depends strongly on the diophantin properties of the classical pe- 
riods T c i v T c u. If the fraction T c i /T c i is commensurate (in this case the classical 
Hamiltonian flow is T c /-periodic) we can describe accurately the behaviour of the 
dynamics on a classical period [0, T c {\ (see theorem 3.12). In opposite, if the frac- 
tion F c \ I T c i 2 is a bad approximation by rationals (we suppose T c ;. / T c ; 2 is Roth 
number) the autocorrelation function collapse in the set ]0, T s ] where T s is order of 
l/h s (see theorem 3.24). For large time we use the continuous fraction expansion 
of T c i t /T c i 2 to analyze some possible periods for linear approximation &\(t) (see 
theorem 3.35). 

In the last section, we use the quadradic approximation a.z(t) (see definition 
4.6) of the autocorrelation function, valid up on a time scale [0, 1/h.P] where /S > 
1. In this quadradic approximation appear three revivals periods T reVl , T reV2 and 
Trev-n or order 1 / h depending on the Hessian matrix of the function F at the point 
(El, £2)- If we suppose T r e Vl ,T r ev 2 and T re v 12 are commensurate, we can proove 
and analyze the revivals phenomenon (see theorem 4.16 and corollary 4.17). In the 
last subsection we compute the modulus of the revival coefficients (see theorem 
4.19). 



2. General points 



2.1. Some basic facts on semi-classical analysis. To explain quickly the philoso- 
phy of semi-classical analysis, starts by an example : for a real number E > 0; the 
equation 

(where A ? denotes the Laplace-Beltrami operator on a Riemaniann manifold (M, g)) 
admits the eigenvectors cpk as solution if 

Hence if h — > + then A/ c — > +oo. So there exists a correspondence between the 
semi-classical limit (h — > + ) and large eigenvalues. 

The asympotic of large eigenvalues for the Laplace-Beltrami operator A» on a 
Riemaniann manifold (M,g), or more generally for a pseudo-differential opera- 
tor P;„ is linked to a symplectic geometry : the phase space geometry. This is 
the same phenomenon between quantum mechanics (spectrum, operators alge- 
bra) and classical mechanics (length of periodic geodesies, symplectic geometry). 
For more details see for example the survey [Labi]. 

2.2. Quantum dynamics and autocorrelation function. For a quantum Hamil- 
tonian P/, : D (P;,) C T~L — > H, where H is a Hilbert space, the Schrodinger dy- 
namics is governed by the Schrodinger equation : 

= WW- 

With the functional calculus, we can reformulate this equation with the unitary 
group U(t) = je~'£ p '< j . Indeed, for a initial state ipo E H, the evolution is 
given by : 

ip(t) = U(t)ip e H. 

We now introduce a simple tool to understand the behaviour of the vector ip(t) : 
this tool is a quantum analog of the Poincare return function : 

Definition. The quantum return functions of the operator P/, and for an initial 
state ipo is defined by : 

r(t) := (ip(t),ipo) n ; 
and the autocorrelation function is defined by : 

a(f) := |r(t)| = |(^(0/*0>«l • 

The previous function measures the return on the initial state ipQ. This function 
is the overlap of the time dependent quantum state ip(t) with the initial state ipo. 
Since the initial state ipo is normalized, the autocorrelation function takes values in 
the compact set [0,1]. 

2.3. The Hamiltonian of our model. For our study, the quantum Hamiltonian is 
the operator : 

P h :=F(P 1 ,P 2 ) 

where F is a polynomial of 1R [X, Y] which does not depend on the paramater h; 
Pi and P2 are the Weyl-quantization of the classical one dimensional harmonic 
oscillator : 

Pj (xi,£i,x 2 ,£>) = &>/ (xj + £j^J II 
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withcv\,L02 > 0. It is well know that the Hermitte functions (e„ /tn ) n := (e„ ® e m ) w meN i 
is a Hilbert basis of the space L 2 (]R 2 ). Let us consider for all integers (n,m) the 
eigenvalues of P\ and P2 : 

T n := wi/z + , p m := co 2 h (m + ; 

so, we get immediatly that for all integers (n,m) 

F (Pi, Pi) ( e n <8 e m ) = F (t„, p m ) (e„ <g> e m ) . 



2.4. The autocorrelation function rewritten in a eigenbasis. Now, for a initial 

vector i/>o = Yj a n,m^n,m we have for alH > 

n,meN 2 

\n,meIN 2 / n,meN 2 



so, for all f > we obtain 



n=0 m=0 



and 



-j-t*J -(-uu 

*(t)= E E Kml 2 *-^^) 

a(f) 



I 2 p-'i f (Wm) 



+00 +00 

E D l fl «,m| 2 e~''' 

n=0m=0 

The aim of this paper is to study this sum, but unfortunately this function is too 
difficult to be understood immediatly. 



2.5. Strategy to study the autocorrelation function. The strategy for simplify the 
sum function t 1— >■ a(f), performed by the physicists ([Av-Pe], [LAS], [Robil], 
[Robi2], [BKP], [Bl-Ko]) is the following : 

(1) we define a initial vector tpo = E a n,m e n,m localized near some regular 

n,meN 2 

Liouville torus of energies (E\, E2) : consequently the sequence (a n ,m) n mgisj 2 
is localized close to a pair of quantum numbers uq, uiq (depends on h and 
on the Liouville torus (E\, E2). 

(2) Next, the idea is to expand by a Taylor formula's the eigenvalues F (r n , ]i m ) 
around the Liouville torus (E\, E2) : 



F(r n ,n m ) = 

3F 3F 

F (r„ , }im ) + hcv t (n - n ) — (t„ , }i mo ) + hcv 2 (m - m ) — (t„ , }i mo ) 

1 3 2 F 1 3 2 F 

+-u\h 2 (n - n ) 2 ^ ( T n ,Hm ) + ^cv z h 2 (m - m ) 2 ^ ( T n ,}hn ) 

+<xi 1 W2h 1 (n - no) (m - m ) (t„ , }i n , ) + ■■■ 

(here r no , }im is the closest pair of eigenvalue to the pair E\, £2)- As a con- 
sequence we get for all f > 



a(f) 



\&n,m\ 

n,»ieN 2 



2 -it \coi(n-n )^(r„ ,fi mQ )+-+coiw 2 h(n~n )(m-m )^ ? (T„ ,ft,„ ) + 



1 



(3) And, for small values of f, the first approximation of the autocorrelation 
function a(f ) is the function 



a x (f 



i2 - 

\,m e 



-if[wi(n-n )^(T„ ^ mo )+a; 2 (m-mo)|f (t„ ,/(„ !() )] 



E \ a "> 

n,iti6N 2 

and for larger values of f , the order 2-approximation is given by 

a2 ^) . = £^ | flnw |2 e -!'4 w i( n - n o)^(Tn ^m )+---+wia) 2 ft(n-no)(m-mo)^(Tn ^ mo )] 
n,mGlN 2 

Insection3, we study in details the function t \— > H\(t) and t i— >• a2(£) in section 4. 
2.6. Choice of an initial state. Let us define an initial vector i/>o = E a n,m e n,m 

n,m£N 2 

localized near a regular Liouville torus of energies E := {E\, £2) where Ei 6 [0, 1] 
and £ 2 G [0,1]. 

Definition 2.1. Let us consider the quantum integers «o = n o(h,Ei) and ?«o = 
m (h, E 2 ) defined by 

n : = argmin|T„ - Ei| ; m := argmin \}i m - £ 2 | . 

n m 

Remark 2.2. Without loss of generality, we may suppose that the integers hq and 
niQ are unique. 

The integer n$ (resp. mo) is the eigenvalues index of the operator from the 
family Pi (resp. P2) the closest to the real number E\ (resp. £2). Since the spectral 
gap of Pi (resp. P2) is equal to co\h (resp. co 2 h ) we have, for h — > : n ~ 

Now, we can give definition of our initial state : 

Definition 2.3. Let us consider the sequence {a. n ,m) n m ez 2 = ( a n,m(h)) nm£Z 2 de- 
fined by : 



#n,m : — ^hX 



h d 2 



n-riQ m-m 
K hX [ ^1 ,"2 



h d 2- 



where the function x is non null, non-negative and belong ot the space S ( K 2 ) . The 
parameters (6'^, 5' 2 ) g]0, 1 [ 2 . We also denote 



h 



X 



Tn Tjiq Mm F^o 



/z<> 2 



£ 2 (N 2 



Let us detail this choice : 
(1) the term x 



localize around the torus (£1, E2) (for techni- 



cal reason we localize around the closest eigenvalues to (E\, £2). 
(2) Constants S[ and S' 2 are coefficients for dilate the function x (the reason to 
take < 6j < 1 is the following : it is the unique way to have a non-trivial 

localization (not tend to {0}) and a localization larger the spectral hi 2> h). 

So, clearly the sequence (a n ,m) n m £ ^ 2 (Z 2 ). Now, let us evaluate the constant of 
normalization Kj,} start by the : 

Lemma 2.4. For a function cp £ <S(1R 2 ) and (£1,62) G ]0,1] 2 then we have uniformly 
for («i,»2) G 1R 2 

1 + Mi s + W2 ' 



E 

2 ,K+"l|>2,|s + n 2 |>2 



£2 



0(e? + e?). 



Proof. We see easily that, uniformly for (u\, u^) 6 K 2 we have 

(I + U\ S + U-2 \ 



e,s&2 2 , \e+ui\>z, |s+«2|>j 



£l £ 2 



0(1). 



Next 



£,s€Z 2 , |€+«i|>2/ l s +"2|>2 



< E 

e,sez 2 , \£+u%\>2, |s+u 2 |> 2 



(^) 



2N 



2N 



t-\-U\ S + «2\ 

£l ' £2 / 

£ + Ml S + U 2 \ 



<P 



And, similary we have 



e,seZ 2 , \£+ui\>i, |s+k 2 |> 2 



ei e 2 

f l + U\ S + «2 
V £l ' £ 2 

i + Ui s + u 2 \ 



r 2N 



(^ + Ml) 2N 



£1 £ 2 



/ 1 + «i S + M 2 \ 



£l £ 2 



To conclude the proof, we apply that to the functions ip(x,y) := x 2N (p(x,y) and 
ip{x,y):=y 2N q>{x,y). □ 

An obvious consequence of this lemma is the following result : 

Proposition 2.5. We get 



1 



+ 0(«°°); 



^3(x 2 )(0,0)h— 
hence ||fln,m||^( N 2) = 1 + 0(h°°). 

Proof. By the Poisson formula and the lemma above we get the equality : 



E x 2 U 



n, mi 



n-n m-m \ L ^ + ^_ 2 



y[-i h s'-i- 

I ,-s 



z(x z )(o,o)+ E 

f,s6Z 2 , K| + |s|>l 

•J'+S'-Zr, 



h 6[-\ h S' 2 -\ 



u> 2 



0J\ 



0)2 



ft^-^ir (o,o) + o(fe~)- 



Now, with the basic equality 

2 / n — hq m — niQ 



E 

n,m6N 2 
-1 +oo 

E E * 2 

n — — oo m= — co 



h S[-l ' h 8> 2 -\ 



2 ( n — Hq m — mo 



and with the lemma above we see easily that 

~ 1 +0 ° 2 f 11 ~ n m ~ m 



2 ( n — uq m — niQ 

2_( 2j x 

n=0 m=—co 



E E x 2 

n=—oo m=— oo 



6 



hoo -1 



2 / n-n m-m 
n=0 m=-oo \ « 1 « 2 



Finally we get 



X 



hence 



n — no m — mq 



( 2 {N 2 



3[X Z ) (0,0)h s 'i +s 2- 2 + O(h 



+ 0(/z°°) 



For finish, we write 



\ a n,m\\e 2 (K 2 ) — K h E 



X w r 



n — hq m — mo 



/A 



Sf* 2 (0,0) + OQT) =l + 0(h°°). 



□ 



2.7. Technical interlude : the set A. In this subsection, we introduce the set A C 
N 2 , this set is useful for making approximation for autocorrelation function. Start 
by the definition : 

Definition 2.6. Let us define the set of integers A = A(h, E\, E 2 ) by : 

A := |(n,m) £ N 2 ; \t„ - t„ | < co^" 1 and \ji m - p mQ \ < a^ 2 } 

= Un,m) £ N 2 ; \n - n \ < h 5 ^ 1 and |m - m \ < h 6 ^ 1 } 
where < 5 t < 1; and we define the set T = T(h, E\, £2) by : 

T := N 2 - A. 
We have the following usefull lemma : 
Lemma 2.7. If we suppose for all i £ {1, 2}, S\ > 5j then we have 



\&n,n 



0(h°°). 



n,m£T 



Proof. The starting point is the following inequality : 



f,m&T n,meZ 2 ,\n-n \>h^ 



\a n ,m\ + E \a n ,m\ 

n,meX 2 , |m-m |>fc' ) 2- 1 



Since the function x 2 is in the space <S(]R 2 ), for all integer N > 1 we have 



E 

n,ra£Z 



n-n 



h 6 i 



-1 



2N 



i n,m I 



m — mo 



2N 



l n,m I 



O(l). 



n,ffl6Z 2 

Without loss generality, we may suppose that Uq = mo = 0. Next we write 

2N 

,2 h 2N(S[-l) 



n,m£Z 2 , |n|>^l _1 



n,m£Z. 2 , |n|>Jt*l _1 



S'-l 



1 

2N 



= o (h 2N ^-^) 
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In a similar way, we get 



because 5- > Sj, this implies V~] |fln, m | 2 = 0(/z°°), so we prove the lemma. □ 

«,mer 



3. Order 1 approximation : classical periods 

3.1. Introduction. In this section, we use a Taylor's formula to expand the phase 
term tF (T n ,}i m ) /h in the variables (n,m) in linear order. In this approximation 
appear two periods T c ; and T c ; 2 of order 0(1) depending on the gradient of the 
function F at the point {E\, E2). 

3.2. Linear approximation and classical periods. 



Assumption 3.1. Here, we suppose that |^ (E\, E2) 7^ 0, |y (^1/ ^2) 7^ 0. 
3.2.1. Semi-classical and classical periods. 

Definition 3.2. We define semi-classical periods r sc / and T sc / 2 by : 



dF 



In 



l sch ■ dF 



and T c 



2/r 



H (T„ , f( mo ) Wl " C ' 2 |f ( T »0' Fm ) ^2 

So, in linear order approximation, we have : 

Proposition 3.3. Let a a real number such that a > 1 — 2mm5j. Then, uniformly for 
allte [0,h a ]: 



+ o (V +2min<5i_1 



r(t) = e- itF ( T '^ l "-o) /h £ | fl „ /m | 2 e 

n,m£N 2 

Proof. Let us introduce the difference e(f) := e(t,h) defined by 



e(t) : = 



-2l7Tf 



E |fln, m | 2 ^~^ F(T "' Fm) - e -' fF ( T »o^-o)/ ft J] |a„, m | 2 e 

n,niEN 2 n.raeN 2 

For all integers (n, m) G N 2 the Taylor-Lagrange's formula (at order 2) around 
(t Mo , ^ mo ) on the function F gives the existence of a real number 9 = 8 (n, m, n§, m§) £ 
] 0,1 [such that 

271/2 (n — no) 27t/j (m — wiq) 



E (^71/ — ^ ("*"i?0' Fmo) 



^^F(p n m ) 2 2 2 , 1 9 2 f (P«,m) 2,2/ x2 



2 ax 2 



a 2 F , 

+ axay < -'°"'" ! ' ) Wia;2/l ( M ~ M °) ( m ~ m °) ' 

withp„, m = p(n,m,n ,m ,h) := (T ng +6(T n - T no ),u mo + 6(u m - u mg )) . 
So, we get 



e (0 



E K 



, 2 e 



-lint 



-ilntK nM {K) 



where we have used the notation 



_ tof(n - w ) 2 8 2 F (p n>m ) - m ) 2 3 2 F (p n>w ) 

t^n,m\") ■ — ; =tttt r 



4tt 



ax 2 



47T 



3Y 2 



+ 



ha>icv 2 {n - n )(m - m ) d 2 F (p U/ „ 



In dXdY 
With the sets T,A and by triangular inequality, we obtain for all f > 



e(t)< 



-2int[ ^31+231 



n, me A 



-i2ntR„ im {h) _ ^ 



' 2 \&n,n 
n,m£T 



For all t > 0, for /z small enough and for all integers (n,m) £ A, we observe that 
^(w-« ) 2 3 2 F(p„ /m ) < tK ^-l. 



47T 



ax 2 



f/zo; 2 (m-m ) 2 9 2 F(p„, m ) < ^^-1. 



47T 



thwiw 2 {n - n )(m - m ) d 2 F (p n ,m) 



In dXdY 
where K\,K2,Ki2 > are constants which does not depend on h. Indeed : let us 
denotes by B ( (E\, E2), r) the Euclidian ball of dimension 2 with center (£1, E2) and 
radius r; since lim/,_>, (T„ ,fi mo ) = (E l7 E2) we obtain that Ve > 0, 3}iq > 0, such 
that for all h < ho, (T nQ ,}i mo ) £ B ((Ei,E2),e); next for all integers (n, m) £ A, we 
have |0(t„ — t„ )| = hw\Q \n — n \ < co\h 3x and \6(}i m — F«o)l = hcv 2 9 \m — mo| < 
(x>2h Sl , this means that for h small enough (h < ho) we have 

Pn,m £ B ((E 1 ,E 2 ),£); 

therefore we obtain for all h < Hq, 



ax 2 



(Pn,n 



< sup 

(jf,y)eB((Ei,E 2 ),e 



^F 

ax 2 



and this quantity is > and does not depend on h. Next we have for all t £ [0, /i"] 

t \Rn,m(h)\ < Kth^ 23 ^ 1 + K 2 h a+23 ^ 1 + K 1/2 h a+ ^ +s ^ 1 

< Mh*- 1 (h 23 ^ + h 232 + h 3 ^+ 3 ^ = ZMh 2 ™* 3 '^- 1 ; 

where M := max (K\, X2, X12) ; with (by hypothesis) 2mini5/ + a — 1 > 0. This 
implies that for all t £ [0, h a ] and for all integers (n, m) £ A we get 

e -i2ntR n ,m(h) _ ^ _ q ^2mirK5 ; +a-l^ . 
and consequently we have for all t £ [0, h a ] 



E 

n.mgA 



&n,m I £ 



-i2ntR n ,m(h) _ ^ 



= j Jj 2 m ini,+«-lj £ | flwm |2 =0 j jl 2mW i+ «-l 

Finally, for all f £ [0, h a ] we have e(t) = O (^min^+a-i) 



□ 



The semi-classical periods T sc /. depend on the parameter h. Later we consider 
two cases : T sc u / T sc i 2 £ Q or not. Consequently we don't make commensurability 
hypothesis on the number / T sc ;, valid up for all h > 0, so we prefer introduce 
two quantities which does not depend on h to make latter commensurability hy- 
pothesis. So we replace semi-classical periods T sc /. by semi-classical periods T c \. 
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Definition 3.4. We define classical periods 7^ and T c / 2 by : 

27T , „ 271 



and T d := 



If (El/ |f (El/ £2)^2 
An obvious remark is that for all / £ {1, 2} we have lining T sc i = T c \.. 

Proposition 3.5. Let rbea real number such that r > — min <S,. Then, uniformly for all 
te [0,h T \: 



E 

n,m6]N 2 



2 



n-iin , m—m ( 



n,ni6N 2 



Proof. We observe that 



\&n,n 

n,m6N 2 



~l~2 \an,m\ 
n,m£A 



'- + ■ 

H ' scl 2 



-2int V^ + V 



< £^ 2|fl„ /H 
n,ra£r 



27Tf(n-n ) [ J---L 



1 1 

27Tf(m - m ) ( - ^ 



because |e ,x ie ! ' x 2 - e ! ' Y V T 2| < 2 |X X - Y t \ + 2 \X 2 - Y 2 \ 
Next for all f > we have 



Intin — no) 



1 1 



2nt{n — no) 



\ 1 sc/jTc/j 



and we know that 



_ 2n H {T„ ,u mo ) - |f (£i, £2) _ 

11 If (El/E2)|f (T„ ,Fm ) 

first, applying the inequality of Lagrange we obtain : 

dF 1 \ dF iv r ^ 

^ {T„ ,U mo ) - — {t lr t 2 ) 



< sup 

*,yeB((Ei,E 2 ),l) 



(§1 



K 2 



< MJ (t„ - Ea) 2 + { } i mo - E 2 f < Mh 



(T„„,f( mo ) - (El, E2) [|r2 



V2 



where M > and does not depend on h. 

On the other hand, since we suppose |f (Ei, £2) 7^ 0, there exists £1 > and 
r\ > such that for all (x,y) 6 B ((Ej, E 2 ), ri) we get 

8F 



ax 



>£l- 



We have seen that hence that there exists hi > such that for all h 6 ]0, /zi [ 

^ mo ) e B((E 1 ,E 2 ),r 1 ); 
as a consequence the application h <— > dF ^ , r is bounded on the open 

■Sx\ E l' E 2)Sx{T:n ,}hn ) 

set ]0, hi [; indeed for all h £ ]0, h\[ 



5X (El, £2) gx ( T no'Fmo3 
10 



1 

< < +00 

e l 



hence, with M' := ^M/z^i 
Next, since 



, for all h G ]0,hi[ we have \T ch - T sch \ < hM'. 



< 



LUILL>2 



dF dF 



< 



sup 

,x,y£B((£ 1 ,E 2 ),l) 



dF_ 

dX 



0/3/) 



< 00 



there exists a constant C\ > which does not depend on h such that for all 
h 6 ]0,/zi[ we get 1/T sc j — 1/T C ; | < C\h. In a similar way there exists C2 > 
and hi > such that for all h G ]0, J^f we g et 1 1 / T sc i 2 — 1 / T c / 2 1 < C?h. As a conse- 
quence, for all h G ]0,/z*[ where /z* := min/i,-, for all f G [0, h T ] with t G M, and for 
all integers (n,m) G A we have : 



t(n-n ) 



< C x h v+S \ 



t(m — mo) 



T sc l 2 T d2 



< C 2 h v+h ; 



we thus obtain for all t, (n, m) G [0,li"j x A 
1 1 



tin - n ) 
Therefore 



t[m — mo) 



T sc i 2 T cll 



< Mh T+mhlSi . 



\ttn,m\ 

n,meN 2 



< ATiMh T+mmd < £ |a„, w | 2 + O (ft 00 ) = O (ft T + mm ^j . 

n,ffl£A 



□ 



3.2.2. Comparison between classical periods and the time scale [0,h K \. In proposition 
3.3 the hypothesis on a is that a > 1 — 2min^,; therefore with S{ G]|,l[ we can 
make a "good choice" for a; i.e. to have a < 0. Hence for /z small enough we 
obtain : 

[0,T d; ] c [0,/z ff ]. 
Next, since — min<5, — (1 — 2min<5 ( ) = — 1 + mint!), < we get 

^— mirnJ, ^ Jjl— 2mirnJ ; . 

this means that we can choose to take T = a. 



3.2.3. X/ze Zz'near approximation aj. In conclusion, the linear approximation of the 
autocorrelation function on the time scale [0, /z a ] is : 

Definition 3.6. The linear approximation of the autocorrelation function is 



ai :f |fln, 

II,fflEN 2 



m I e 
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3.3. Geometrical interpretation of classical periods. The periods T c \. have geo- 
metrical interpretation. For E\,E2 > consider the energy level set Me u e 2 := 

p^ 1 (E\) n P2 1 (E2) C K 4 , this manifold is isomorphic to the torus y^^S 1 x 

tt^S 1 , here S 1 is the one-dimension circle. Start with the calculus of the Handl- 
er 



tonian flow of p = E(p\, p 2 ) with an initial point niQ £ Me v e 2 - So the Hamilton's 
equations are 







( \ 


x 2 (t) 










—ax\(t) 


V 6(0 J 




K -MO / 



where we have used the notation a :— 4§ (E\, E2) co\, b := |y (El/ E 2 ) w 2- For all 
/ £ {1,2}, let us consider the complex number Zj(t) := Xj(t) + i£j(t); from the 
Hamilton equations we obtain the equalities Z\ (t) = —iaZi(t), Z 2 (f) = — z'£>Z 2 (f). 
Therefore we get 

Z 1 (0 = Z 1 (0K lflt , Z 2 (f) = Z 2 (0)e- iht 

and 



2£ 



|Z!(0)| 2 = x?(0) + g?(0) 
this means that the Hamiltonian's flow in complex coordinate is given by 



fl I,|Z 2 (0)| 2 = x 2 2 (0)+^(0) = ^; 



1 y 



Zi(0) 
Z 2 (0) 

In angular coordinate the flow is given by 

01,0 \ / #1,0 

02,0 J V 02,0 

argZ^O) 
2tt 

ian's flow : 

27T 27T m 271 

1 



Z a (f) 

z 2 (0 



f fl 



In 



with 0y o — 



1] . So we have exactly the classical periods of the Hamilton- 

2tt 



|| (Ei,E 2 ) w\ 



ff (El/ E 2 ) <^2 



It's well know that if the periods are commensurate the flow is periodic on the 
torus. In opposite the flow is quasi-periodic on the torus. 

3.4. The principal part of the function aj. Now, let us study in details the func- 
tion a\(t) on the time scale [0, max T c i] . Start by a technical proposition : 

Proposition 3.7. For allt>0 we have 

y~l \&n,m\ 



-2mf| ^ + ^2 

c/i do 



n,m£2 



Six 2 ) (0,0) 



a'! 



f - 



f 



h*2 



U'2 



l d 2 - 



Proof. The trick here is just to use the Poisson formula, so let us consider the func- 
tion Of defined by 

R 2 -> C 



Of : < 



-V2 I 
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where t £ IR is a parameter. For all integers (n,m) £ Z 2 we have 



, -2int \ -t— t T 
|«n,m| e v c 'i c ' 2 y = n f (n,w). 

So clearly, the function Qt £ iS(R 2 ), then the Fourier transform $ (Of) is equal, for 

all £i,£ 2 eR 2 

$ (O f ) (ft, ft) = / / n f (x a , x 2 ) e - 2 ™&e- 2mx & d Xl dx 2 ; 



CO J —CO 



therefore for all ft, ft G IR 2 we g et 



e -2i7t(wo?i+mo&) f 2 \ ( ( t \ h 3 !" 1 ( „ t 



It comes from the Poisson formula the equality 

£ n t {n,m)= £ y(n,)(*,s) 

which gives the proposition. □ 

Since the function # (^ 2 ) £ <S(]R 2 ), we observe that only index I, s £ Z 2 such 
that I + or s + ^y- are close to zero are important in the sum. More precisely : 

Definition 3.8. For all t > 0, let us define the integers £j(t) = £i(t,h,E) as the 
closest integers to the real numbers — f/T c ;.; i.e : 

e i (t) + ^-=d(t,T cl .Z); 
where d{.,.) denote the Euclidiean distance on ]R. 

Remark 3.9. Without loss of generality, we may suppose the integers ii(t) are 
unique. On the other hand, for all integer I £ Z such that i ^ £j(t) we get : 



Tcli 



1 

> -. 
- 2 



Lemma 3.10. Uniformly for t > we have : 
Proof. Since $ (x 2 ) £ 5(R 2 ) we have 

six 1 ) (ft, ft) 



\/k,d £ N* 2 , 3B k4 > 0, Vft,ft e R 2 , 



< D /c,rf 



(i + |ft|) fc (i + lftlf 



Next, it then follow from the proposition above and from the lemma 2.4 that for 
all t > 



-1 
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Next, for all £ > 



a i(0 = E \ a w 



,2 + 



n—nn , m—rfir 



E 

n,m6Z 2 -N 2 



_2i7tt 53+5=31 

\ u n,m\ e 



thus 



ff(x 2 )(o,o) 



-d(T ch Z,t), — — d(T ch Z,t) 



L0 2 



For finish, we observe 



< E k,m| 2 + o(/z°°). 

«,meZ 2 -N 2 

| Hn,m | 

n,m6Z 2 -N 2 



+ co -1 -1-1 -1 +oo 



n=0 m=-oo 



n — — oo m— — oo 



n=-oo m= Q 



and an obvious consequence of the lemma 2.4 is that meZ 2 -N 2 l a "-" ! I = ^ (^°°) ■ 

□ 

3.5. Behaviour of the function ai : case G Q- In this subsection we suppose 



T^ = l eQ; hence aT ch = bT clr 

Definition 3.11. If the classical periods T c i ,T c u are commensurate the classical 
period of the global system is defined by T c ; := aT c u = bT c i . 
Now, we can formulate an important result of this section : 

Theorem 3.12. We have : 

(i) for f real such that f G T d Z we get (i.e. for all i G {1,2}, d (f, T d .Z) = 0) 

ai(0=l- 

(ii) If there exists i G {1, 2} such that d (T c/ .Z, f) > fr 1 ^/ then : 

a t (0 = O(n- 

Proof. The first point (i) is clear. For the second : it follows from the lemma 3.10 
that 



ai(£) 



1 



W)(o,o) 



d (T dl Z, f) , d (T ch Z, t)\+0 (h 00 ) ; 



L0 2 



since the function £ (x 2 ) G 5 (1R 2 ) we have 



Vr? G N, 3D, > 0, V^fc G M z , 3 (* Z J 
and therefore 



< 



Ij^-l I,<5 2 -1 

-d(T ch Z,t),-——d(T ch Z,t) 



U> 2 



< 



Da 
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Thus, if there exists i 6 {1/2} such that d (T C /.Z, f) > cofi 1 s i then there exists 
£ > such that d (T C ;.Z, f) > cOjh 1 ~ s 'i~ £ and thus for all ^ G N we obtain 



hence we get 

5 fx 2 



-d(T ch Z,t), —d(T ch Z.,t) 



U>2 



-d(T ch Z,t), —d(T ch Z,t) 



a> 2 



0{h° 



□ 



3.6. Behaviour of the function ai : case ^ £ Q. Let us now tackle an important 



case : the case 



— is not a fraction of Q. Here, there not exists classical 



common period, the Hamiltonian flow is not periodic on the torus. 

First, we note that, in view of lemma 3.10, the behaviour of the function ai is 
given by the function : 



1 1 y 



Six 2 ) (ao) 



six 



-d(T ch Z,t), —d(T cl ,Z.,t) 



a>2 



Therefore, since the function £ (x 2 ) belongs to the space S (M 2 ), we need to explain 
simultaneously the evolutions of the distances d f T e ^.Z, depending on time. In 
another formulation, we want to analyze the behaviour of the Euclidian distance 
between the segment line (OMf) where O := (0,0), Mt := (f,f) an d the lattice 

Tc/jZ x T C ; 2 Z depending on the time t and the number ^ = |. Precisely, we want 

to compare the distance d ((OMf) , T c / Z x T C ; 2 Z) with the real number h & i 1 . For 
example, if for a time t* the distance is larger than h i we get ai (f*) = 0(h°°). 

Start this new subsection by some geometrical results and latter we explain the 
study of the autocorrelation function ai(t). 

3.6.1. Some general points. In angular coordinates the Hamiltonian flow is : 

[0,1] 2 -> [0,1] 2 



32,0 



Without loss of generality we may suppose that the initial data is ^ ^ 1/0 ^ 



q I and that a,b <0. Therefore the Hamiltonian flow is given by q>t 



af 
bt 



where we have used the notation a := — j£- > and b := — ^ > 0. Recall that 



2tt 



function 



and T, 



I 2/r I 

m 



So, to understand the behaviour of the 



s(f)(ao 



-d(T c/l Z,f),--— d (T d2 Z,t) 



we need to explain the evolution of the Euclidian distance d (<pt,Z 2 ) depending 



on time f and on the real number j . 
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3.6.2. Suppose \ verify diophantin condition. J. Liouville proved in 1884 the follow- 
ing theorem : 

Theorem 3.13. (Liouville). For all algebraic irrational number 6 with degree d > 2 
there exists a constant C = C(6) > such that the inequality 

C 



> 



holds for all rationals £ . 

H 

In other words, algebraic numbers are bad approximation by rationals. Finally 
in 1955 K. F. Roth has considerably improved this result (he was awarded the Field 
medal in 1958). 

Definition 3.14. We say an irrational number 8 satisfy a e-diophantin condition 
(e > 0) if and only if there exists a constant C E > such that 



> 



holds for all (p,q) 6 Z x N*. We denote by C E the set of irrationals 8 that holds 
e-diophantin condition. We say that an 6 irrational 8 is a Roth number if and only 
if 8 e f]C e ;le 

£>0 



Ve > 0, 3Q > 0; V(p, q) G Z x N*; 
There is a lot of Roth numbers examples : 



> 



r,2+£- 



Theorem 3.15. (Thue-Siegel-Roth). Every real algebraic irrational number of de- 
gree d > 2 is a Roth number. 

We have also the (see for example [Cas]) : 

Theorem 3.16. The Lebesgue measure of Roth's numbers is infinite. 



Remark 3.17. Let 8 a £-diophantin number. Since for all p we have 1 8 - 
\0 - p\ < \ < ^r; thus we obtain that < Q < \ holds for all e > 



> C e and 



r2 ._ 



Now, we estimate the Euclidian distance between the set {(pt} te \Q 7] and Z^ 
Z 2 -{(0,0)}. 

Notation 3.18. Let us denote by A and T the following orthogonal lines 

A := Vect(aei + be 2 ), T := Vect(-bei + ae 2 ) 

where (ei,e2) is the canonical basis of the vector space K 2 . Let us also considers 
7Ta the orthogonal projector on the line A and ny the orthogonal projector on the 



line T. 



Here we suppose 8 = ^ is a Roth number. 



Lemma 3.19. For all e > there exists < K £ < C E , here C E denotes the Roth constant 
of 8 = §, such that 



> 



\\n A (ne 1 +me 2 )\\ R 2 
holds for all (n,m) £ Z 2 . 



\nei + me 2 \\^2 E ' 



\n T (nci + me 2 )\\ ]R 2 > 



\ne x + me 2 \\ 1 E 2 £ 
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71 



=5 ( b ) ^ e u ™ tar y vector or the 



Proof. Let us denotes by u = ( U ^ 
line A, so we have 

||7T A («ei + me 2 )\\ R2 = \{u,ne 1 + me 2 ) B i\ 
= \nui + mu 2 \ = |«i| \n + m9\ ; 
since 6 is a Roth number, for all £ > there exist C £ > such that 

||tt a {ne x + me 2 )\\ K 2 > \u t \ t^t+i 



> 



\u x \C e 



\tiet + me 2 \\^z £ 



In a similar way we ge 



Q 



\n T (ne! + me 2 )\\ R 2 > \u 2 \ n > 



I "2 1 C £ 



n| 1+£ " \\ne 1+ me 2 \\]+ £ ' 



□ 



therefore with K E : = min ( | U\ \ C e , \ U 2 \ C e ) < C £ we obtain the lemma. 
A consequence of this lemma is : 

Theorem 3.20. For all e > there K e < C e , here C £ denotes the Roth constant of 
9 = §, such that for all f > 



l+£' 



Proof. We observe that for all f > the point q>t belongs to the line A, thus there 
exists a pair (nt, tnt) £ Z 2 such that 



d(<p f ,Z 2 ) = Chpt ~ {n t e\ +m t e 2 ) 



R- 



> \\n T [n t e 1 +m t e 2 )\\ R2 , 
applying the lemma above we get 



\n t ei + m t e 2 \\^ 2 £ 



On the other hand we have the majorization 

Oft ~ (wf^i + wi t e 2 
and, by triangular inequality we obtain 

\\n t ei +m t e 2 \\ E 2 < 

— ± 

Therefore, since 



0<p t 



R2 - 2 



R2 2 



Ocp t = t Va 2 + b 2 , we get for all i > 0, £ > 

R 2 ° 



d (<pt,Z.l) > 



1+6 " 



□ 



Corollary 3.21. For all e > and for every i-j s 
rf f<P(,Z 2N | <« f > ==1 — 



0, 



x/2 1+E 
2 


c 


0,^ 




i 


V2 


V 7 > 




2 
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Proof. Suppose d (ft(0),Z*) < r\, it then follows from the theorem above that for 



all e > there exists a constant K E £ 



0,i 



K r 



holds for all t > 0; i.e. 



such that 

—TH <r l 



: — + tVa 2 + b 2 . 



□ 



Remark 3.22. Since r/ 6 



0, 



c 



0,^ 



we have ( -f 



1+5 > \/2 



Notation 3.23. For e > and t] > 0, let us denote : 



1 



Va 2 + b 2 



'/ 



^ 



with r/ small enough such 



Theorem 3.24. For all e > 0, for every r\ £ 

that f,j (e) > max T c /. and for all k > 1; there exists a constant > which does 
not depend on /z such that the inequality 

|ai(f)| < D k h k( - 1 - max ^t]- k 
holds for all t £ [max T c ; ., f J? (e)]. 

Proof. Our starting point is that for all f > max T c ;. we have 

d(f,T c ,.Z) = d(t,T cI .N*). 



Next, since 



I — I 



and T, 



c/ 2 



I27TI 



we get 



d(t,T ch Z) = d(af,N*) 7 d(t,T ch Z) = d(bt,N*). 
Therefore, from the corrolary above (by contraposed) we obtain for all f £ [max T c ; , f ^ 



d ((a*,bf),i 



> 



and since the norms || (x,y) || R 2 and \x\ + \y\ are equivalent on R 2 , there exists a 
constant C > such that 

d(f,r c/l Z)+d(f,T ci2 Z) > Crj 

holds for all t £ [max T c j., f )? (e)]. 

Next, since the function # (^ 2 ) belongs to the space <S> (K 2 ) we have 



Vfc£N 2 , 3M/t>0, V£i,£ 2 £]R 2 , 



thus we obtain for all t £ [max T c ;., (e)] 



5 (a: 2 ) 



< 



(ICil + lfcl 



[x 



hi' 2 -l 

d(t,T dl Z),-——d (t,T ch , 



U>2 



< 



M k 



"'1 



-d (t,T ch Z) 



C02 



-d(t,T ch Z] 
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max(cOi) k Mk 



maxfafM* 1 =ma , .^M^Mi-max^)) -k 



□ 



Applying this theorem with t] = h s where the real number s belongs to ] 0, 1 — max [, 
we have : 

Corollary 3.25. For all e > 0, s 6 ] 0, 1 — max 8', [ and for h small enough such that 
ih s ( £ ) > max T c i', the following equality 

ai(0 = O(n 

holds uniformly for all t £ [max T c / , ^s(e)]. 

Notes on time scales. From a pratical point of view, we must verify that for all 

e > 

t h s(e)<h« 

where a > 1 — 2 min (£,•) . Indeed we have : 

Proposition 3.26. Suppose min ^ > |,/or all e > and for all s £ ] 0, 1 — max <$( [ we 

W " 2Va 2 + b 2 V 
Proof. For all s > and for all £ > , since 

*M-7!»TS> (<*>**"* -t) 

for rz — > we have the equivalence 

f/is ~ D £ h~~ 

, i 

where D £ := , _ _ (K E ) 1+£ > 0. On the other hand, we see that for all £ > we 

V 3 z +b z v 

have 1 — 2min 8j < —j^. — ■ Hence 

h^r- <fti-2nun<y i> 



Fherefore, we obtain 



1 / 1 max^-l 

1 | (K f )TTi^- 2mm ^ - — I < 1 [ I^l-zminJ,- 

2 +b 2 y ' 2 ) - v^+v v 2 



V2' 



□ 



Use of continued fractions. Now, we can wonder what are the accurate times 
when d (cpf,Z*) < r\ ? To solve this problem we will use the continued fraction 
theory. 
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Some useful theorems. The continued fractions are essentially used for the approx- 
imation of real numbers. There exists two types of continued fractions: the finite 
continued fractions representing rational numbers and the infinite continued frac- 
tions representing irrational numbers. For all irrational number 6, there exists a 
pair sequence (q„, p n ) eN 2 such that 

El < — 
q n ~ ql 

holds for all n > 0. This sequence is given by the continued fractions algorithm 
(see [Ro-Sz], [Khi]). Geometrically speaking, the construction principle for this 
sequence is as follows (see [Arn2]) : consider v := (0,1) and v_\ := (1,0). It is 
obvious that these points lie on different sides of the line y = 6x. By induction : let 
the vectors zty-i and v k be constructed whereas to construct the new vector c>/ c+ i, 
we add to the vector z^-i the vector v k as many times as we can in such a way the 
new vector z^+i lies on the same side of the line y = 6x as the vector : 

Vk+l = a k v k + v k _i 

i.e 

qk+i = a kqk + qk-i 

Pk+i = a k pk + Pk-i 
where (ak)k>o * s a se q uence of integers strictly > 0. 

We note that the sequence (q n ) n is strictly increasing. With the standard nota- 
tion continued fractions we have 

[flO/01/- ■•/«*/•••] := «0 + - — ; — —[ — 



«1 + - T i 

«2 + ^ 



end we have (see for example [Khi]) the relation [uq, a\, . . . , a„ s 

Example 3.27. The number n is given by : n = [3, 7, 15, 1, 292, 1 . . .] . 

Approach time. Let us denote by D the line y = jX = 0x. Hence, for all f > we 
have : 

[Oft] C D. 

For a fixed n > we wish to find the point M n of D such that d (M„, {cj n , p„)) = 
d (D, (q n , p n ) ) . In other words, we wish to find the time T n such that d (f T „, {qn, Pn)) 
d (D, (q„,p„)) . 

Proposition 3.28. For all n > 1 the unique r n > such that d (f Tn , (q n ,p n )) = 
d(D, (q n ,p„)) is given by 

_ aq n + bp„ 
Tn ~ a 2 + b 2 • 

Moreover we have 

d (f Tn ,Zl) < * — < — . 
v ' Va 2 + b 2 qn q n 

Proof. For n ^ fixed, we want to find t > such that d (f t , (q n , p n )) = d (D, (q n , p n 
This means that we want to find t > such that 



Oft -L ( (q n ei + p n e 2 ) -Oft 
i.e. : to find t > such that 



Oft, {q n e\ + p n ez) -Oft ) = 0. 
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Consequently we solve the equation at(q n — at) + bt(p n — bf ) = and we find for 
non-null solution : f = a ^^;" . Therefore, at time t = t„ := a ^'^ b ;" we obtain 



{<pT n ,{<tn,Pn)) = d 
1 



a 2 q n + abp n abq„ + b 2 p„ 
a 2 + b 2 ' a 2 + b 2 



, (q n , Pn) 



b 2 



b 2 (ap„ - bq„f + a 2 (bq„ - ap„f. 



Since for all integer n we know that 
all integer n, so we deduce that 



</</ 



< 4^-, i.e. \q n b — p„a\ < ^ holds for 



'iti 



d (cp Tn , (q n ,p n )) < 



1 



a 2 ,a 2 



'b 2 



2 ^ < 



1 1 

< 



a 2 +b 2 y q\ q 2 Va 2 + b 2 qn qn' 
For conclude, we note that 

d(cp Tn ,(q n ,p n )) > d {q>T n ,^f 

holds for all integer n. 



□ 



We wish to generalize this result : we wish to analyze the behaviour of the 
distance between the set Z 2 and the flow q>t when Ms in a neighbourhood of the 
time t„. 

Notation 3.29. For r > let us denotes by B(r n , r) the closed ball of center r n and 
radius r > : 



B(r n ,r) : = 
Proposition 3.30. For all r > 



aq n + bp n _ aq n +bp n 
a 2 + b 2 ' a 2 + b 2 



2 'a* 



d (* Z ^ a^F V ( O + " (a2 + b2) J + l^ +rb(a2 + ^ 
holds for all t € £>(T„,r). 

Proof. We begin with the following inequality : for all t > and for all n 

d (<p f ,Z 2 ) < d (<p f , (q n ,Pn)) ■ 
Next, it's clear that for all £ 6 B{x n , r) we have 



a 2 q n + abp„ \ / abq„ + b 2 p„ 



Therefore, for all t £ B(t„, r) we have 

d{(ft,{q n , Pn)) 



V a 2 + b 2 



,rb 



< 



a 2 ^„ + abp„ 
a 2 + b 2 



ra-q n 



( abq n + b 2 p„ 
{ a 2 + b 2 



rb - p n 



^ V (abpn - + ra(a 2 + b 2 )) + (.b* - a 2 P „ + rb(a 2 + b 2 ))^ 
1 



< 



a 2 + b 2 



2 'a 2 



? ) ( — + 



□ 



21 



Remark 3.31. For r = we obtain 



d[cp tr Zi) < 



b 2 



< 



Va 2 + b 2 <?n 



and we obtain again the result of the proposition 3.28. 

Now, let us give an asymptotic equivalent (for n — > oo) of the real number T n : 



Proposition 3.32. For n — > oo we have r n ~ Oq,,; rw/tere Q 
Proo/. We just write the fraction r n / Clq n : 

x n _ aq n +bp n a 2 +b 2 a 



a+be 



bp,, 



> 0. 



Qq» a 2 + b 2 q n (a + b6) a + b6 q n (a + b6) 
and since lim^^oo p n /q n =6 we obtain that limj^oo T n /dq n — 1. 



□ 



Now, let us come back to the autocorrelation function approximation a\. Start 
by a notation and a remark : 

Notation 3.33. For ji > let us denotes by Af, = A], (0,ji) the following set : 



A h := [q n G N, ^„ G 



^miniJj— 1— j< ^jl— 2mirn5,+/( 



C 



^miniJJ— 1 ^1— 2min<J,' 



Renwrfc 3.34. Fory > 0, we have of course 

If we suppose that the set Ah is non empty, we have some periods for the func 
tion ai, indeed we have : 



Theorem 3.35. Suppose A), 7^ 0, then 

sup |ai(T„) 

n€{meN,(? m (e)GA} 



O(^). 



Proof. Applying the Taylor-Lagrange formula on the function [x, y) 1— > $ (x 2 ) (x, y) 
near the origin : for all f > there exists 6 = 6 (t, h, T c u, T c u) 6 ]0, 1[ such that 




d(T ch Z,t) 



Li>2 



-d(T ch x,t) )=d(x 2 ) (0,0) 



3ff (* 2 ) 














L ^1 



d(T ch Z,t) t 6 d(T ch X,t) 

CO 2 



h s' 2 -i 

d(T dl Z,t),9 d(T ch X,t 

C02 



We know that for all n > 1 the distance between the part of the flow cp Tn and the 
set Z 2 is strictly lower than Hence, if we suppose that from a certain point, 

like n > N , the sequence 



is strictly lower than h s (with s > 0), then we 
obtain the majorization d (<p Tj! ,Z 2 ) < h s . Therefore, for all i 6 {1,2} we have 
also d (r„, Trf.Z) < h s . Consequently for all i 6 {1,2} we get h 6 i~ 1 d (t„, T c ;.Z) < 
h s 'i- 1+s . Since we suppose s > — min S', + 1 + y with j,t > we deduce that : 

tf f i- l d{T n ,T d ?) <tt l ; 
and, for ft small enough, we obtain 



03\ OX 



h s[-\ h s> 2 -i 

-d (T dl Z, Xn) ,6 d (T d2 Z, ttn) 



C02 
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d (T c/2 Z,t„) 



-d(T dl Z,T»),0 



d(T c/2 Z,tT„) 



< Nft?*; 



where M, M' > are constant which does not depend on ft. Next, it comes from 
the Taylor formula written above that 

sup |ai(r„) 



II 



sup 

ne{meN,q m (6)eA h } 



h«L- 



L0 2 



-d (T ch Z.,T n 



-d (T dl Z, t„ 



dx 



h s' 2 -i 

d (r £ j Z,T„) ,6 d (T d Z,t„) 



3t/ 



d (T d Z, t„) , d (T C ; 2 Z, T„J 



□ 



Counting of the sequence q n . In view of the theorem above, let us now tackle an 
important problem : what is the cardinality of the set Ah ? Note that since the 
sequence (q n (8) ) ngN is strictly increasing we have 

#{A h } =#{nen, q n (6)EA h }. 

Start by a simple majorization of the integer # {Ah} '■ 

#{4}<#{Nfl [H^-r-^^s^ <E[S(h)] + l 

where 5(h) : = fei-Zmin^ _ ^min^-l-p and denotes the integer part of x. 
Then, for ft — > we have the equivalence <S(ft) ~ ftl-2jnini j +^ < g Q we g et a ma _ 
jorization of the integer # {^} in order }i 1 - lmmS i+V , 

Nevertheless, find a minoration of the integer # {.4;,} is more difficult; but it's 
cleat that for all n* > 1 there exists ft* 6 ]0, 1[ such that : 

Vmm^-l-, U *l-2mm^+^ p ^,,(0)^ = {^(0)}. 

In order to estimate the integer # {.4;,} , we must know the distribution of the se- 



^min^'— 1— ft jjl—2mmt>i+}i 



quence (q n (8) ) neN on the real axis (in particular on the compact set 
depending on the number 8. Let's try to give some distribution examples of the se 
quance (<7„(0)) neN . 

An exemple : the golden ratio. The golden ratio cp is the unique real roots of X 2 - 
X — 1 = 0, i.e. f = (1 + y/5) /2. The continued fraction of the golden ration is : 

1 



) 



1,1 1,... 



1 + 



Consequently the golden ratio is that one of the most difficult real number to ap- 
proximate with rationals numbers. An another particularity of the golden ratio is 
that the sequence of the denominators (q„) „ from the continued fraction algorithm 
is equal to the Fibonacci sequence (T n ) n : 



We note that 




F 



n+1 



lim 

H-S-+00 F„ 
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Next for n — > +00 we have also : 




We have also the following property : 

Proposition 3.36. Denote by q n (x) the sequence of denominators from the continued 
fraction algorithm of the number x;for all 8 6 R, n > we have 

q n (6) > F„. 

In the general case for any 9 irrational number, we have the following theorem 
(see for example [Khi]) : 

Theorem 3.37. (Khintchine-Levy, 1952). Almost surely for 9 G 1R we have 

lim q n (6)n =K; 

n->+oo 



where K denotes the Khintchine-Levy constant K := e 121n P) > 1. 
Thus for instance from a certain point we obtain : 

\kJ < qn {9)< (Ik 

The study of the distribution of the geometrical sequences ^2^) ) M anc ^ ( ( 1^0 " 

on the compact set / l mm5 ;- 1 -/ 1 / / ? 1 ^ 2mm ' 5 /+/' i s easy; unfortunately it does not 

provide accurate informations on the distribution of the sequence (qn(9)) neN ■ 

Open question. Do we know the denominator distribution of (q n (0))neN with the 
real axis depending on 9? More specifically, for a non-empty compact set of diam- 
eter 5 > includes in is it possible to estimate the number of elements of the 
sequence (q n {9)) n£N in this compact set depending on the numbers 9 and S ? 

4. Second order approximation : revival periods 

4.1. Introduction. Our next aim is to use a more accurate approximation of the 
function 1 1— > a(f). In this section, we use the quadradic approximation 82(f) of the 
autocorrelation function, valid up on a time scale [0, l//z^] where f> > 1. This ap- 
proximation is a consequence of a Taylor formula on the the term tF (t„, ji m ) /h in 
order 2. In this quadradic approximation appear three revivals periods T reVl , T reV2 
and T reVl2 (of order 1/h). 

-\ 2 p ~\2 p 

Assumption 4.1. In this section, we suppose ^? (Ej,E2) 7^ 0, ^ x5y (Ei,E2) 7^ 
0, g(Ei,E 2 ) ^0. 

4.2. Quadradic approximation andc revival periods. 

4.2.1. Semi-classical revival and revival periods. 

Definition 4.2. Let us define the semi-classical revival periods T sreVl , T sreV2 and 

Tsrevi2 by • 

T - 4/T _ 47T 

S '' roi h PL (T ~ )0J 2' ls '^ '- h £F (T W 2 ; 

" gx 2 v "0' rm ) w \ n gy2 I '"c rm ) w i 

T 4/r 

1 srev 12 — d2F " • 

24 



So we get the approximation : 

Proposition 4.3. Let f> a real number such that f> > 1 — 3mmSj. Then we have uni- 
formly for t G [0, /z p ] : 

e +itF { T »oV»'a) /h r(t) 

"0 (»-» ) 2 ("--»'q) 2 ("-»o)('«-'"q) \ 

>2 Z'srePj Tsrei'2 T^srevyi J 



y~! \ttn,m\ £ 
n,meN 2 



-2l7Tt 



"-"0 _|_ m-mn 



+0 ftj0+ 3min *«— 1 



Proof. The principle is the same as in the proof of proposition 3.3. Here we use 
the Taylor-Lagrange formula at order 3 : for all pair in,m) G N 2 there exists 9 = 
9 (n, m, no, mo) G ]0, 1 [ such that 

F (t„, u m ) = F (r no ,u mo ) + d I^^co 1 h (n - n ) + dF a> 2 h (m - m ) 



1 9 2 F(T„ ,^ mo ) , 2 



+7: 



2 

1 a 2 F 



3X 2 



07j/z (n — no) 



dy 

2 1 d 2 F (Tn ,U mo ) 21 



dY 2 



cvpi (m — mo) 



IdXdY 



(Tn ,U mo ) wxwih 1 (n - n ) (m - mo) + - 



ld 3 F(p nim ) 3 



ax 3 



(vfh (n - no)" 



ld 3 F(p n>m ) 2 ; 3/ \2/ \ , 13 3 F(|0«,m) 2,3/ w \2 

'-a^a^n (n - n ) (m - m ) + - d ^dY 2 lW2 ^ ~ H °' ^ ~ ™ ' 



2 dX 2 dY 



+ - 



1 d 3 F(p„ rtnj 
6 3X 3 



co\h 3 (m - mo) , 



withp n<m = p(n,tn,n Q ,m Q ,h) := (t„ + 0(t„ - T„ ),u mo + 9(u m - }i mo )). 
Next, we observe that for all pair (n, m) G A and for all f G [0, /z p ] 



t (n - n ) 3 h 2 



< fcfH^-l. 



t (n — uq) (m — mo)h* 



2 7 

f (n — no) (m — mo) h 



t(m — mo) h 



< ^p+ZSi+Sz-t 



hence, since /3 > 1 — min (3<!>i,2<5i + $2,<$i + 2^2,3<$2) = 1 — 3min<5,', for all t G 
[0, ft"] and for all pair (n,m) G N 2 we get 

e ~2int((n~n fh 2 +(n~n ) 2 (m~m )h 2 +(n-n )(m-m ) 2 h 2 + (m~m ) 3 h 2 ) 



= 1 + ^/jP-l+3nun< 5 i 
And the statement of the proposition is established. 

For the same reason as in definition 3.4 we introduce the revival periods : 
Definition 4.4. Let us defines the revival periods T reVl , T reV2 and T reVn by : 

An „ An 



□ 



h^ I {E 1 ,E 2 )cv 2 ' ^ h^(E lr E 2 )cv 2 2 
An 



1 rev n 



h WZ5Y ( E l' £ 2)^1^2 



Clearly for all j G {1,2,12} we have lim^o Fsrevj/Frevj = 1. The three semi- 
classical periods T srev . depend on h as well as their quotients. Since we will con- 
sider period quotients afterwards, is it preferably to study revival periods than 
semi-classical revival periods;for that we use indeed : 
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Proposition 4.5. Let v a real number such that v > —2 min 5j. Then we have uniformly 
forte [0,h v \: 



\a-n,n 



-lint 



, -2int( ""-'"n) 2 + <"-"Q»"oA 

n-n m-m (»-»p) 2 . (»i-m ) 2 ("-"pX"--"^) \ 

'sdj 'sd 2 ' J 'rrai ' J 're»2 Tr '' I 'l2 / + Q f /l" 



+min(<5] / (5 2 ) 



Proof. The principle is the same as in the proof of proposition 3.5. With the parti- 
tion N 2 = A II T and by triangular inequality we have 



\ttn,m\ 

n,mEN 2 



-linti "-"0 | "'-'"0 I '"r"0 )2 I ("'-'"0> 2 I '"-»o)( m - m O) '\ 
\ ' J sc/ 1 T sd 2 ' srev l ,srev 2 T srev K J 



2i7it( n ~ n ° I I ( "'" o)2 I ( "'' m " )2 I '"-"o)' m -"'o) Y 

^ T scf^ T scl 2 Tm '\ Trl!v 2 T rev\1 J 



+2 \a n>n 

n,m£A 



2nt(n - n Y 



< 2 \a n ,n 

n,m£T 
1 1 



Tsrevi 



2nt(n — no)(m — mo) 



2nt(m — ntoY 
1 



Tsrev2 Trev2 



1 srev\i 



reviz 



because \e iX - e iY \ < 2 \X - Y\ 
Next we observe that 



Tsrevi Trevi 



47T 

to\h 



fj£ { T n ,Um ) fp {El,E 2 ) 



First we have 



^2 (,tl,t2 J - \ T n ,}lm ) 



< sup 

(x,y)£B((E 1/ £ 2 ) / l) 



V (g) (*,y) 



R 2 



< Ma/ (Ei - t„ ) 2 + (E 2 - f( mo ) 2 < Mfr 



(El/ £2) - (T«o,^mo)llR2 



V2 



2 ' 



where M > is a constant which does not depend on h. 

On the other hand, since we suppose (Ej, E2) 7^ 0; there exists £\ > and 
ri > such that for all(x,y) G B ((Ei, E 2 ), ri) we get 

(*/y) > £ i; 

and we have seen that there exists h\ > such that for all h £ ]0, /;i [ we have 

u mo ) £ B((E 1 ,E 2 ),r 1 ); 

therefore the application 

1 



h 1 ^ 



fx? (El, E 2 ) |^ (t„ , }l mo ) 
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is bounded on the open set ]0, h\ [, indeed for all h £ ]0, h\ [ we have 

1 

< -j < +00; 



(El/ £2) ( T n /^m ) 



"1 



and with M' := 2tt^ f ora llfc £ \%h\\ we obtain |T srWl - T ret>1 | < M'. 



Next, since 



<h 2 



§x^ (El/ E2) f^l ( T « / /-'mo) 



167T 2 



< Khr 



, there exists a constant Q > 



where K := ^ sup (ie/y)6B({ElA)/1) | fg(x,y) 
(which does not depend on h ) such that for all h £ ] 0, h\ [ we have 1 1 / T sr£ . ni — 1 / T reni | < 
Q/z 2 (eg. take C\ := KM). In a similary way : there exists C 2 , Q 2 > such that for 
all /z 6 ]0,fr 2 [weget |l/T sre » 2 - 1/T re „ 2 | < C 2 /i 2 and for all h £ ]0,h u [ we get also 

\l/T S ren 12 ~ 1/T r en n \ < Cl2^ 2 - 

Next, for all t > 0, for all pair (n,m) £ A and for all /z < min (/*Zi,fo 2 ,?Zi 2 ) we 
have 

< Qlilh 2 ^; 

< Ci 2 |t|/J ,5l+ ' 52 ; 



f (n - n ) 

t (« — Mo) (m — 7«o) 

f (m — moY 



Tsrevi Trevi 



Tsrev2 Trev2 



< C 2 \t\h 2S i; 



hence for all f £ [0, h v ] where v is a real number such that v > —2 min Sj we obtain 



t(n-n ) 2 C_J_- J_ 
t (n - no) (m - m ) 



t(m — mo) 
1 



1 



1 



Tsrev2 



Tsrev^ T r evi2 



rev2 

< Mh v+2minSi 



where M := 3 max (Q, C 2 , Q 2 ) . 

So we proove that : there exists a constant M > such that for all h < min (h\, hi, h\2) 
and for all v > —2 mini?, we get 



y~! \d-n,m\ 
n,m£N 2 



n-n ) 2 (m-m ) 2 ("-"q)('"-'"q) \ 
Tsrev-y Tsretfn Tsrev^j J 



2»7rff"~"° 1 1 ( "'" o)2 1 ( "'~ m « )2 1 ( "-"o) (m -"'o) Y 

\ T sc/ 1 T sd 2 Tret, l T re v 12 J 



<2 £ \a„, m \ 2 + Mh v+2mklS ' £ \a n>m \ 2 . 

n,m€:T n,m£A 

It follows from the lemma 2.5 (the lemma 2.5 says En,mer l a «,m| = O and 
from 

E l fl 'V"| 2 < \ a «,m\ 2 = 1 + O (h°°) 

n,me& n,meN 2 



that 



I ^n, m I 

n,meN 2 



-2l7lf 



k-kq m-m (H r «p) 2 . (»'-m ) 2 ("-" )(m- mo) \ 
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Q ( ^u+2min<5; 



□ 

4.2.2. Comparison between revivals periods an the time scale [0, h & \. Recall here that 

the parameters (s' if £] \, 1 [ 2 with J- > Sf, recall also that the real coefficients «,jS 

would require a > 1 — 2min^, and f> > 1 — 3minJ;. Next we observe that 1 — 
3min^; — (1 — 2min<5;) = — min^, < and — 2min<5 ( - — (1 — 3min Sj) = min^- — 
1 < 0; hence for h small enough we obtain 

j^l— 2min<J, ^ ^1— 3trun<5; ^ j^— 2min<$, 

So we can make a "good choice" for parameters a, f> and v : indeed we can choose 
a, f> and v such that v < < — 1 < a < 0, therefore (for h small enough) we have 

[0,T di ] C [0,h K ] c [0 r T reVi ] C [0,^] C [0,h v ). 

4.2.3. The quadradic approximation a.2- So, the quadradic approximation of the au- 
tocorrelation function on the time scale [0, h.P] is : 

Definition 4.6. The quadradic approximation of the autocorrelation function is 



4.3. Revival theorems. 



2 



t Yh-»- m-mp (»-„p) 2 (m-mp) 2 ("-"pXm-Mq) \ 
V T sdj T scl 2 Trev l Trev 2 Trev 12 J 



4.3.1. Preliminaries. 
Resonance hypothesis. 

Definition 4.7. We say that the revival periods T reVl , T reV2 , T reVl2 are in resonance 
if an only if there exists ( — , — , — ] 6Q 3 such that 

Pl _P2 T _ Pl2 T 

qi qi qn 

Notation 4.8. In this case, we introduce the notation Tf rac := ^T reVl = ^T reV2 = 
^T reVl2 . And for all ;' G {1,2} let us also consider the numbers ry := p\ 2 qj, s j '■ = 
qi2Pj , and clearly for all £ {1,2} we have T reVj = j:T reVl2 . 

Preliminaries. To make progress in our study we need to introduce a new func- 
tion \p c i with two artificial variables t\, t 2 . 

Definition 4.9. Let us define the pseudo-classical function ip c i : 

m—mQ 



lpcl(h,t 2 )--= E \ a n,m\ e sd l 



2 



n,meN 2 

So we get the obvious following property; first the function xp c i is doubly-periodic 

(i) for all pair t lr t 2 > we have f d (fj + T sdl , t 2 ) = ip c l (h, t 2 ) ; 

(ii) and for all pair t\, t 2 > we have also ip c i (t\, t 2 + T sc ; 2 ) = tp c i (t\, t 2 ) . 
This function have no immediate physical significance, but if the time f i and t 2 

are equal : 

(iii) for all t > we have \p c i(t, t) = ai(f). 
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Some lemmas. 

Notation 4.10. Let us consider the sequence {Q n ,m)n,m = { Q n,m{v\,q\,P2,qi,PYi,qYi,h)) ni 
with (n,m) G Z 2 defined by : 

0„ m e - 2i7r (fr("-»o) 2 +|g(»-"o)('n-m )+ffi(m-mo) 2 )_ 

The periodicity of this sequence is caracterised by the following easy proposi- 
tion. 

Proposition 4.11. For all p\, q\, P2, qi, P\2, q\2 G Z, f/ze sequence (6 n ,m) n m verify 

^n+£i,m = ®n,m 

if 'and only if the integers £\ and li satisfy the following equations : 

V(», m) EZ 2 > + 2 J^n + S [1] 

V(n, ») £ Z 2 > + ^2 m + ^2 n s „ 

Example 4.12. An obvious solution is £i = t/iSi and £2 = q2 s 2- 

For two periods i\, £2 G Z 2 let us consider the set of sequences l\, ^—periodic 
with his natural scalar product. 

Definition 4.13. For a fixed pair li,l z e (Z*) 2 we define 6^ 2 (Z) the set of se- 
quences t\, I2— periodic in the following sense : 

&i v i 2 (Z) := lu n ,m G C z2 ; Vn, m G Z 2 , u n+ ^ m = u n ,m and u, hm+ e 2 = u n , m \ ■ 
So we have the elementary : 
Proposition 4.14. The application 

(U,v) H> {U,V) 6 := Yl E u n,mV n ,m ■ 

n=0 m=0 

is a Hermitean product on the space &£ ^ (Z) . 
We have also the obvious following remark : 

J- „ 2inkn 2mpm 

Proposition 4.15. Let us consider <p n 'f m := e 1 e c 2 where (k, p) G Z 2 ; then the 

family < ((pnm) ( zs an orthonormal basis of the space vector 

[\ ' /n,meZ?) fc=o..^ 1 -l,p=0..^ 2 -l 

4.3.2. The main theorem. In the following theorem we show that the function t t-> 
82(f) near the period T^ r(JC can be written as a finite sum of ip c i with arguments 
shifted. Indeed we have : 

Theorem 4.16. Suppose resonance hypothesis holds; then there exists a family of 
£\ + £-2 complex numbers (depends on h) : ( c fc 1/ fc 2 ))c 1 6{o £j_i}/c 2 6{o 11. where the 
integers £\, £2 G Z 2 are solutions of equations from proposition 4.11; such that 

a2 ( f + I/ rac J = E !C C h,*2$cl * + Tfrac + Y Tsc h' t + T frac + J- T scl 2 ) 
7 fc 1= fc 2 =0 V <1 *2 / 
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+0 ( }i a + 2minS i- 1 



holds for all t 6 [0, h a ]. The numbers c^fa are called fractionnals coefficients; and 
for all ki G {0...£i - 1}, k 2 G {0...£ 2 ~ 1} 



with = b kl/h (h) = (a h ,<p' 



Proof. Let us denote the integers n 
function e(t) : 



n — hq, m := m — mo and consider the 



e(0 : = 



a 2 U + l/ rflc J - E c kik$cl ( * + 7/ rflc + Y^scly 1 + 7/ rilc + ^-T sc ; ; 



fc 1= 0fc 2 =0 



n,mGN 2 



, 2 e' 
m 1 <- 



-2int^— -ZinT, 



frac T 



-2int~, 



-2inT 



1 e 



:i 2 e 



frac 7 } 



-2intJ- -2in^- -2int^- -Zin^i -lint 



v i e ''1 e 



'2 e 



12 e 



her 



~2in'- 



' E E C h,k z ^cl [ t + Tfrac + Y Tc h' f + T f"ic 
k 1= 0k 2 =0 \ 1 2 



i2 e 

rue ~r e i c l 



n,m6N 2 



2int^— -2iTzT fracr H- -lint^- -2mT fmcl ^- 
g J sdj g J sd 2 g J sc( 2 



-ZintJ^ ZintJ^ Zint^a- 



e ™i e rco 2 e 

"EE Chfalpcl [ t + Tfrac + yT sdl , t + Tf rac + yT sch 
/ Cl =0fc 2 =0 \ 11 ll 

Since the sequence {0 n , m ) n m G &£.£ 2 (Z) with £j = qiS\ and ^2 = ^2 there exists 
a unique decomposition of the sequence (6 n ,m)„ m on the basis | (Vn'jw 
; indeed we have : 

0?z,m = E E ^ > k^,k 2 ( Pn,rn' ■ 
k l= 0k 2 =0 



jt=0..i 1 -l,p=0..i 2 -l 



where 



0, ^l'^ \ Therefore we get 



E(t) 



ZiTtt-r*- -2inT. 



t x -\l 2 -\ 

E E E \ a n,m\ Kk 2 e Tsdl e 

n,meN 2 k 1 =0k 2 =0 



-Zint-n 



-ZinT, 



: 'i e 



:l 2 g 



frac j 



-2intj£ 2int^ 2171(7422- k fc 

g ^g w, g i^^KiAa 

' E E C k t ,k 2 tycl ( f + T/rac + y^sdy f + T/rac + Y^sch 



, -2(^(7^- -2inT fmcl J!- -2int^a- -2inT fracT at- 

h, , o sdj g - sc'i g sc( 2 g J J sd 2 



E E E ^iA-2 e 



n^iGlN 2 fci=0A: 2 =0 



-2i7rtT^ 2t7t^ 2171(7^- k fc , 
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E 3 

k 1 =0k 2 =0 n,m£]N 2 



-2(7Tf T ^ 2inT fracl ^- - 2in k x f- 

E E C ^ 2 E e le ie 71 



-7-iut-^— -2iTzT fmcl M- _ 2ink * 

sd 2 g J scl 9 „ 



■'2 e 



lint^fi— -2inT, 



4-14-1 

E E E l«n,«r*W Tsd ie 

n,mEN 2 fc a =0fc 2 =0 



free t. 



-2int^— ~2inT 



He 



:< 2 g 



frac' 



g Jreoj g Jre» 2 g 'reojj g tj g e 2 
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a 1 ■ -2inki-P- -2i7tk 2 f- 1 r 2 — j j ,t , 

And since Cfc a/ fc 2 e 1 e 2 = %Jfc2 e e we deduce tnat 



£ (0 



1-14-1 
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11,1I16N 2 fc 1= 0fc2=0 
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E E E \ a n,mWxXi e SC ' le SC ' le SC ' 2e SC ' 2e Te 



-lintjt ZintJ^ 2;'7rf T =- 
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To finish, we use the partition N 2 = A II T and we just consider indices in the set 
A for the sum. Hence there exists constants C\, C 2 , C\ 2 > which does not depend 
on h, such that 



(n - n ) 2 



<-rev 1 



(n - n )(m - m ) 



(m - m ) 2 



T, 



rev 2 



holds for all pair (n, m) £ A and for all t £ [0, h a \. 

As a consequence for all t £ [0, h a ] and for all pair (n, m) £ A we get : 



-2intJl 2intj£ Tint-^- 



"i e 



Trcih 



2 6 



12 - 1 = O [h 



a+2miri(J, — 1 



□ 



If we take t = we obtain : 
Corollary 4.17. Under the same hypothesis we have 
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4.4. Explicit values of modulus for revival coefficients. Our final aim is to com- 
pute the modulus of revival coefficients. The idea is to split the sum | in two 
simple parts. These parts look like that Gauss sums, but in fact with a little dif- 
ference. Here we propose a simple way to compute this sums and we don't use 
sophisticated theory. Start with a notation and a remark : 

Notation 4.18. For I > 1 and for integers p et q such that p A q = 1, let us consider 
for all integer k E {0...£ — 1} the following sum 



' n=0 



E~2in^m 2 2Ankm 
e i e~r- 



Therefore for all k G {0..1 - 1} 

14(4 p,<?) | = 



Theorem 4.19. Suppose resonance hypothesis holds and suppose also ^ G Z; 
then we obtain : 

|^!,/c 2 | 2 = |4i (h,PlS\Jl)\ 2 |4 2 (^2/P2S2/4)| 2 - 

Proof. For all jfcj G {0...^ - 1}, k 2 G {0..i 2 ~ 1} we have 
1 



1 

4^ 



E 

(n,m)eZ,/e 1 Z,x 
2ink 2 m 

m=0 



e V '/i '/12 12 / e 'i e '2 



?2 > e fi e '1 V 1 112 / 



E* 

n=0 



and since l\ = = q\q\iP\ we obtain 



1 



%,k 2 



1 2inknm p<? 5 ^1 1 - . 0-1 ? linn n \ 

\ e ' 2 e 12 > e ?i e n 

m=0 «=0 



For y G {1, 2}, let us consider ^ the following characters : 

Z/^-Z C* 



-linj- 
a i-s- e / ; 



as a consequence we get 
,2 1 r- 



n*2 (x,y)eZ/iiZxZ/t 2 Z 



Xi ( pisix 2 - x(fc x - puqmy) X2 ( P2S23/ 2 - % 



E 

(z,t)ez/l{Zx. 



Xi ( -P1S1Z 2 + z(kx - puqiPit) X2 ( -Pisit 2 + k 2 t 



p. p. 



X\ (pisi [x l -z l \ -ki(x-z) + pu<jiPi(xy-zt) 



X2 (P2S2 (y 2 - f 2 ) -k 2 (y-t) 



Now, because ^ G Z then £j = qiqnPilqiPiPu hence for all x,y,z,t we have 
PnqiPi( x y — zt) G i\Z . Therefore for all x, y, z, t we have Xi (,PnqiPi{ x y ~ zt)) = 
1. Hence 

\K,h\ = 
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f2e2 ^ 
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To finish we can compute 

|4i {h,V\ s \>h)f |<4 2 (h, P2S2, h)\ 2 ■ 
with the following results (see for example [Lab2]) : 

Proposition 4.20. For all pair p, q with p A q = 1 and q odd, then for all k G {0. 

we get : 



•<?-!} 



I P/<?)l 



And 



Proposition 4.21. For all pair p, q with p A q = 1 and q even, then for all k £ {0. 

zwe get : 
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(3 2 

z'/ ^ is odd f/ten |djt(<?/ p, q) \ 



^ if k is even 

else; 
// k is pair 



- else. 
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